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The Dynamical Dipole Mode in Fusion Reactions with Exotic Nuclear Beams 
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We report the properties of the prompt dipole radiation, produced via a coUective bremsstrahlung 
mechanism, in fusion reactions with exotic beams. We show that the gamma yield is sensitive to the 
density dependence of the symmetry energy below/around saturation. Moreover we find that the 
angular distribution of the emitted photons from such fast collective mode can represent a sensitive 
probe of its excitation mechanism and of fusion dynamics in the entrance channel. 
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Production of exotic nuclei has opened the way to ex- 
plore, in laboratory conditions, new aspects of nuclear 
structure and dynamics up to extreme ratios of neutron 
(N) to proton numbers (Z). An important issue addressed 
is the density dependence of symmetry energy term in the 
nuclear Equation of State, of interest also for the proper- 
ties of astrophysical objects [ll, [2, S| ■ By employing heavy 
ions collisions at appropriate beam energy and central- 
ity the isospin dynamics at different densities of nuclear 
matter can be investigated [3, S H S S| ■ 

Here we discuss isospin effects in dissipative collisions 
at low energies, between 5 and 20 MeV/n, were unsta- 
ble ion beams with large asymmetry will be soon avail- 
able. The starting point is that in this energy range, 
for dissipative reactions between nuclei with different 
N/Z ratios, the charge equilibration process in the en- 
trance channel has a collective character resembling a 
large amplitude Giant Dipole Resonance (GDR). Sev- 
eral microscopic transport simulations like semiclassi- 




Time- 
[13 or 



cal Boltzmann-Nordheim-Vlasov (BNV) 
Dependent Hartree-Fock (TDHF) [H Qj 
Constrained Molecular Dynamics (CoMD) [15| have con- 
firmed this scenario. Several experiments based on a 
comparison between two reactions with different entrance 
channel charge asymmetry, leading to the same final 
products, have confirmed the presence of this effect, see 
|16| . [17| and refs. therein. 

The gamma yield resulting from the decay of such 
pre-equilibrium isovector mode can encode information 
about the early stage of the reaction 0, [13, U^ [iM [3 ■ 
This collective response is appearing in the lower density 
intermediate neck region, while the system is still in a 
highly deformed dinuclear configuration. It is therefore 
of interest to look at the influence of density dependence 
of symmetry energy below saturation upon the excita- 
tion and dynamics of the prompt dipole mode. The cor- 
responding emission rates can be evaluated, through a 
"bremsstrahlung" mechanism, in a consistent transport 
approach to the rection dynamics, which can account for 
the whole contribution along the dissipative noncquilib- 
rium path, in fusion or deep- inelastic processes [i9|. In 
this way the data can be directly used to probe the isovec- 



tor part of the in medium effective interaction below sat- 
uration density. 

Here we discuss also the expected emission anisotropy 
of such prompt dipole radiation, due to the neu- 
tron/proton oscillation along a definite symmetry axis. 
This is of interest for planning new experiments and for 
the relation to the lifetime of such transient collective 
mode. 

Exotic Beams 

In the following we shall study the features of the pre- 
equilibrium dipole considering the reaction ^^^Sn +^^ Ni 
("132" system) at lOMeV/A, as referred to the same re- 
action induced by a ^'^'^Sn beam ("124" system). From 
simple arguments a larger initial dipole moment will 
trigger higher amplitude isovector oscillations increasing 
the chance of a clear experimental observation. During 
the reaction dynamics the dipole moment is given by 
D{t) = ^X{t), where A = N + Z, and N = Ni + iVa, 
Z = Zi + Z2, are the total number of participating nu- 
cleons, while X{t) is the distance between the centers of 
mass of protons and neutrons. We note that the initial 
dipole {t = 0: touching configuration) attains a value 
around 45 fm for the exotic ^^^Sn beam, to be compared 
to the smaller value 33/to for the stable "124" system, 
which can be considered as a reference partner in an ex- 
perimental comparison. 

We have employed a mean field transport approach, 
based on the BNV equation, which properly describes 
the selfconsistent couplings between various degrees of 
freedom. The potential part of the symmetry energy, 
Esym/A{pot): 
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is tested by employing two different density parametriza- 
tions, Isovector Equation of State (Iso-EoS) of the mean 
field: i) ^ = 482 - 1638p, {MeVfm^), for "Asysoft" 
EoS, where Esym/A{pot) has a weak density dependence 
close to the saturation, with an almost flat behavior be- 
low po; ii) a constant coefficient, C = 32MeV, for the 
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FIG. 1: Dipole Dynamics at lOAMeV, h = 4/m centrality. 
Exotic "132" system: (a) Time evolution of dipole moment 
D(t) in real space; (c) Dipole phase-space correlation (see 
text). Panels (b) and (d): same as before for the stable "124" 
system. Solid lines correspond to Asysoft EoS, the dashed to 
Asystiff EoS. 



"Asystiff" EoS choice, where the interaction part of the 
symmetry term displays a hnear density dependence. As 
shown in details in refs. [3, 0] these choices represent 
two classes of widely used effective interactions, that still 
require some confirmation from new independent observ- 
ables. The isoscalar section of the EoS is the same in 
both cases, corresponding to a compressibility around 
220Mey. In the numerical simulations a test particle ap- 
proach with 200 gaussian test particles per nucleon has 
been employed. In this way we get a good description of 
the phase space occupation, essential for the low energy 
reaction dynamics. In the collision integral in medium 
nucleon- nucleon (N-N) cross sections are considered, [20| . 
We perform calculations for three impact parameters: b 
= 0,2,4 fm, to cover the region where fusion is mostly 
observed. In order to reduce the numerical noise we run 
twenty events for each set of macroscopic initial condi- 
tions and the displayed quantities are the averages over 
this ensemble. 

In Figure 1 we report some global information concern- 
ing the dipole mode in entrance channel. The time evo- 
lution of the dipole moment D[t) for the "132" system at 
h = Afm is represented in Fig.l (a). A similar behavior is 
seen at 6 = 0, 2 fm,. We notice the large amplitude of the 
first oscillation and the delayed dynamics for the Asys- 
tiff EOS related to the weaker isovector restoring force. 
We can also evaluate the quantity DK{t) — {-^ — -^), 
the canonically conjugate momentum of the X{t) coor- 
dinate. The phase space correlation (spiraling) between 
D{t) and DK{t) is reported in Fig.l (c). It nicely points 
out a collective behavior which initiates very early, with a 
dipole moment close to the touching configuration value 
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FIG. 2: (a) panel: Exotic "132" system. Power spectra of 
the dipole acceleration at 6 = 4fm (in c^ units), (b) panel: 
Corresponding results for the stable "124" system. Solid lines 
correspond to Asysoft EoS, the dashed to Asystiff EoS. 



reported above. This can be explained by the fast forma- 
tion of a well developed neck mean field which sustains 
the collective dipole oscillation in spite of the dinuclear 
configuration with a central zone still at densities below 
the saturation value. 

The role of a large charge asymmetry between the 
two colliding nuclei can be seen from Fig.l,(b,d) panels, 
where we show the analogous dipole phase space trajec- 
tories for the stable ^'^'^Sn +^^ Ni system at the same 
value of impact parameter and energy. A clear reduction 
of the collective behavior is evidenced. 

In order to appreciate if these differences can be ob- 
served experimentally we estimate the gamma yield in a 
bremsstrahlung approach [1^, |2l| : 
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where E^ — fiiu is the photon energy and | D"{lo) \'^ 
is the Fourier transform of the dipole "acceleration" 
D"{oj) = //"°" D"{t)e^'^*dt. For each event to represents 
the onset-time of the collective dipole response (phase- 
space spiraling) and tmax the "damping time", i.e. the 
time step corresponding to an almost flat D(t) behavior. 
In refs. |12l . Il4| the pre-compound dipole radiation 
yield is evaluated via a stationary oscillatory model for a 
time-dependent charged source, described in a TDHF ap- 
proach. We note that the "bremss" procedure used here 
is more realistic in the general case of non-equilibrium re- 
action dynamics: in fact while for stationary conditions, 



when I D"{uj) 



to 



D{uj) p, we recover the same 



results, see [19], in presence of damping effects the tran- 
sient nature of such fast collective mode is consistently 
accounted for ^^. As we will clearly see in the follow- 
ing, the damping dynamics, underestimated in a TDHF 
scheme, is in fact important for the excitation function 
of the fast dipole radiation and for the expected angular 
anisotropy. 

In Fig. 2 (a) we report the power spectrum, | D" {oj) p 
in semicentral "132" reactions, for the different Iso-EoS 
choices. The gamma multiplicity is simply related to it. 



see Eq.Q. We clearly observe a lower value of the cen- 
troid, as well as a reduced total yield, in the Asystiff case, 
due to the weaker restoring force for the dynamical dipole 
in the dilute "neck" region, where the symmetry energy 
is smaller [3,] . Slightly wider distributions are obtained in 
the Asysoft case, due to the larger neutron evaporation, 
that damps the collective oscillation. The corresponding 
results for the stable "124" system are drawn in the (b) 
panel. As expected from the larger initial charge asym- 
metry, the Prompt Dipole Emission is increased for the 
exotic n-rich beam. 

From Eq.(2) we can get the total, energy and impact 
parameter integrated, yield for the two systems and the 
two Iso-EoS. We find 3.0 lO'^ (2.5 lO^^) for ^"^^Sji and 
5.7 10-3(4.4 10-3) for 132^^ in the Asysoft (Asystiff) 
case. 

A detailed analysis of the sensitivity of the results to 
the symmetry energy choice can be performed just fit- 
ting the dipole oscillations by a simple damped oscillator 
model, D{t) = D{to)e'^'^«+''/'''^\ where D{tQ) is the value 
at the onset of the collective dinuclear response, wq the 
frequency, that depends on the symmetry energy choice, 
and T the damping rate, related to two-body N-N colli- 
sions and neutron emission. The power spectrum of the 
dipole acceleration is given by 






(3) 



which from Eq. ([2]) leads to a total yield proportional to 
ujot{ujq + 1/r )-D(io) — ^oTD{to) since ujqt > 1. We 
clearly see the effect of the Iso-EoS on the total yield, 
through the quantity ljo'^t^ that is slightly dependent on 
the system. Hence, from the above relation, the differ- 
ence of the yields associated with two different systems, 
that is the quantity usually exploited in the experimental 
analysis |17|, depends on the Iso-EoS and the sensitivity 
is amplified when using exotic, more asymmetric beams, 
due to the factor D{to)^, allowing for a clear experimental 
observation. It is worthwhile to mention that, according 
to our fit, we find that that the parameter -D(io) rnay be 
less than the touching point dipole amplitude, especially 
in the Asystiff case and for the exotic neutron-rich sys- 
tem. A delay in the onset of the collective response is 
expected and so a more reduced -D(io) with respect to 
the initial "geometrical" value, see also the following. In 
fact, we find that the ratio of the total, impact parameter 
integrated yields obtained with the two Iso-EoS (Asysoft 
relative to Asystiff) is larger in the ^^'^Sn case. We obtain 
1.2 for the ^'^'^Sn reaction and 1.3 in the -^^"^Sn case. This 
result points to other interesting Iso-EoS studies that can 
be performed from an accurate measurement of spectrum 
and yield of the prompt dipole radiation. 

Anisotropy 

Aside the total gamma spectrum the corresponding an- 
gular distribution can be a sensitive probe to explore the 



properties of preequilibrium dipole mode and the early 
stages of fusion dynamics. In fact a clear anisotropy 
vs. the beam axis has been recently observed [23]. For 
a dipole oscillation just along the beam axis we ex- 
pect an angular distribution of the emitted photons like 
W{0) -^ sin^ 6* - 1 + 02^2(0056*) with 02 = -1, where 9 
is the polar angle between the photon direction and the 
beam axis. Such extreme anisotropy will be never ob- 
served since in the collision the prompt dipole axis will 
rotate during the radiative emission. In fact the devia- 
tion from the sin 9 behavior will give a measure of the 
time interval of the fast dipole emission. Just for com- 
parison with statistical compound nucleus GDR radia- 
tion we remind that in the case of a prolate nucleus with 
a collective rotation, for the low energy component we 
can have an anisotropy parameter 02 — —1/4, averaging 
over all possible rotation angles and all possible orienta- 
tions of the collective angular momentum (orthogonal to 
the beam axis), l2J|. Orientation fluctuations can even 
reduce such anisotropy, |2a |. 

These results cannot be translated directly to the case 
of the dynamical dipole. As we see from our calculations 
(Fig. 1) the preequilibrium oscillations extend over the 
first 250 — 300/m/c. During this time interval, depending 
also on the centrality and energy, the deformed nucleus 
may not complete a full rotation on the reaction plane. 
Let us denote by 0^ and (j)f the initial and final angles 
of the symmetry axis (which is also oscillation axis) with 
respect to the beam axis, associated respectively to exci- 
tation and complete damping of the dipole mode. Then 
A0 — (pf — (j)i is the rotation angle during the collective 
oscillations. We can get the angular distribution in this 
case by averaging only over the angle A0 obtaining 



W{e) ~ 1 - (- + -x)P2icos9) 



^sin{(pf-4>i) 



(4) 



where x = cos{(j)f+(j)^) ^y_^^l 

It is easy to see that for A(/) = and A(/) = 27r we 
recover the two cases discussed above. Moreover if 0/ = 
'/'i = 00 (i-e. the orientation is frozen at an angle (pQ with 
respect to the beam axis) Eq.Q gives an 02 = — (1 — 
|sm2(/)o), with a change of sign for (/)q > 55°, i.e. a 
decrease of W{9) around 9 « 7r/2. The point is that 
meanwhile the emission is damped. 

Within the bremsstrahlung approach we can perform 
an accurate evaluation of the prompt dipole angular dis- 
tribution using a weighted form where the time variation 
of the radiation emission probability is accounted for 



W{9)^Y.(3,W{9,^,) 



(5) 



We divide the dipole emission time in At^ intervals 
with the corresponding <^i mean rotation angles and the 
related radiation emission probabilities Pi = P{ti) — 
P(i,_i), where P{t) = j! \ D"{t) p dt/Ptot with Ptot 
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FIG. 3: "132" system, (a) panel: time dependence of the 
rotation angle at b=2fm (dashed line) and b=4fm (solid line), 
(b) panel: time evolution of the emission probability P(i), 
see text, for h = 4fm impact parameter, (c) panel: weighted 
angular distributions for b = 2fm and b = 4fm centralities 
for different symmetry term choices. Dashed lines for the 
Asystiff choice, solid for Asysoft. The Iso-Eos effects on the 
rotation angle are negligible. 



given by P{tmax), total emission probability at the final 
dynamical dipole damped time. 

In Fig. 3, (a) panel, we plot the time dependence of the 
rotation angle, for the "132" system, extracted from dy- 
namical simulations at 6 = 2/to and b = Afm. We note 
that essentially the same curves are obtained with the 
two Iso-EoS choices: the overall rotation is mostly ruled 
by the dominant isoscalar interactions. Symmetry energy 
effects will be induced by the different time evolution of 
emission probabilities. 

This is shown in Fig. 3, (b) panel, for the "132" system 
at 6 = 4/m impact parameter. We clearly see that the 
dominant emission region is the initial one, between 50 
and 150 fm/c, while the dinuclear system rotates of about 
20 degrees, roughly from 40° to 60°. Another interesting 
point is the dependence on the symmetry energy. With a 
weaker symmetry term (Asystiff case) the P{t) is a little 
delayed and presents a smoother behavior. As a conse- 
quence we can expect possible symmetry energy effects 



even on the angular distributions. 

This is shown in the (c) panel of Fig. 3, where we have 
the weighted distributions, Eq.([5]), for b — 2/to and 
b = Afm impact parameters, with the two choices of 
the symmetry energies below saturation. For more cen- 
tral collisions, due to the small rotation of the oscillation 
axis, the delay effect in the asystiff case is not affect- 
ing the angular distribution. For more peripheral reac- 
tions we see a larger contribution at forward/backward 
angles, although the bump around 7r/2 is still present 
due to the decreasing emission probability at later times 
when the larger rotations contribute. Altogether we get 
wider "dipole" angular distributions with respect to the 
beam axis, in agreement with the first available data, 
[231 • Moreover, as evidenced by the results at b = 4 fm, 
we expect to see a sensitivity to the slope of the symme- 
try term below saturation in presence of larger rotation 
velocities, i.e. in fusion events with high spin selection. 

Summarizing we have shown, within a mean field 
transport approach, that in fusion with exotic nuclei an 
enhanced preequilibrium dipole emission can be observed 
with a peculiar angular distribution related to its early 
emission. The features of this collective mode are sensi- 
tive to the density dependence of symmetry energy below 
saturation. The larger emission with the exotic beam will 
enhance the observation of Iso-EoS effects. 

The angular distributions are also sensitive to the fu- 
sion dynamics and dipole excitation mechanism and life- 
time. The results presented here can be important to 
plan new experiments. 

In conclusion the dynamical dipole mode appears to 
be a suitable probe to test the symmetry energy term in 
the nuclear EoS as well as to scrutinize the early entrance 
channel dynamics in dissipative reactions with radioac- 
tive beams. 
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